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ABSTRACT 


This  report  describes  a  system  that  can  provide 
precision  aircraft  navigation  based  on  Loran-C  and  digital 
barometric  altimetry.  The  system  is  expected  to  be  capa¬ 
ble  of  fulfilling  the  critical  accuracy  requirements  of  the 
terminal  area  as  well  as  serving  the  needs  of  area  naviga¬ 
tion.  A  digital  barometric  altimetry  technique  is  de¬ 
scribed,  which  has  sufficient  accuracy  to  provide  the  third 
dimension  for  all  phases  of  aircraft  flight.  Further,  sev¬ 
eral  auxiliary  functions  are  described  that  are  made  prac¬ 
ticable  by  high-accuracy  positioning.  These  are  on-board 
determination  of  the  point  to  release  an  object  for  free 
fall  to  the  ground,  a  built-in  approach  and  departure  sys¬ 
tem,  and  automatic  direction  finding. 
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PREFACE 


This  report  records  the  current  state  of  the  art  in 
the  use  of  Loran-C  for  accurate  aircraft  navigation  and  for 
other  aircraft  functions  that  are  made  practicable  by  accu¬ 
rate  navigation.  During  the  past  five  years  the  state  of  the 
art  has  been  influenced  substantially  by  a  number  of  devel¬ 
opments,  due  princ:  pally  to  the  efforts  of  the  author  and  a 
number  of  his  associates,  namely,  Thomas  A.  McCarty, 
Thomas  W.  Jerardi,  and  Ronald  G.  Roll.  In  addition  to 
their  valuable  contributions  throughout  the  report,  Mr. 
McCarty  has  specialized  in  the  targeting  aspects  of  the 
problem  and  in  the  understanding  of  the  Loran  positioning 
service;  and  Mr.  Jerardi  and  Mr.  Roll  have  contributed 
extensively  to  the  navigation  concepts,  ballistic  solution, 
and  coordinate  conversion.  In  particular,  Mr.  Jerardi 
and  Mr.  Roll  are  responsible  for  Appendix  B,  and  Mr. 

Roll  for  Appendix  C.  Without  their  assistance  and  the 
assistance  of  many  others  during  the  last  five  years,  the 
present  report  could  not  have  been  written.  In  particular, 
we  wish  to  express  our  appreciation  to  Maurice  G.  Vincent, 
William  A.  Mayhew,  Jr. ,  and  Donald  P.  Martini  for  their 
valuable  guidance  and  support.  Publication  of  this  report 
was  sponsored  by  the  Defense  Special  Projects  Group. 
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Navigation  is  perhaps  the  oldest  applied  science, 
i-oug  before  the  Romans,  hardy  men -of  the  sea  were  suc- 
c ''s.M-i'Iy  na  igating  their  ships  to  distant  ports,  and  camel 
.iis  navigated  the  featureless  deserts.  Their  prob- 
however,  were  a  bit  different  than  those  of  a  jet 
vircraft  terminating  its  flight  onto  a  200-foot  wide  runway 
at  160  knots.  To  solve  these  problems  and  those  asso¬ 
ciated  with  traffic  congestion,  navig  tional  aids  are  re¬ 
quired  tnat  permit  the  aircraft  to  determine  its  position 
frequently  and  precisely. 

The  Long  Range  Navigation  known  as  Loran-C  is 
one  of  many  aids  to  navigation,  which  are  now  available 
to  solve  present-day  pro!  lems.  Except  for  a  few.  military 
applications,  very  few  aircraft  use  the  existing  Loran-C 
service.  One  purpose  of  this  report  is  tc  describe  how 
Loran-C  can  be  used  in  aircraft  for  accurate  horizontal 
positioning  —  in  fact,  sufficiently  accurately  to  fulfill  all 
horizontal  requirements  during  landing. 

Since  the  earliest  days  of  aircraft,  altitude  has 
been  estimated  from  measurements  of  atmospheric  static 
pressure.  The  radar  altimeter  was  a  natural  invention, 
but  this  device  has  found  only  limited  acceptance  for  a 
few  special  functions  such  as  verifying  ground  clearance. 
They  lack  the  general  utility  of  a  barometric  altimeter 
wit. -out  which  no  aircraft  flies.  Another  purpose  of  this 
report  is  to  describe  the  high -accuracy,  barometric 
altimetry  techniques  that  have  also  been  applied  in  a  few 
military  systems. 

The  final  purpose  of  this  report  is  to  describe  some 
auxiliary  functions  that  high-accuracy  positioning  makes 
practicable,  such  as  the  on-board  determination  of  the  point 
at  which  to  release  an  object  for  a  free  fall  to  the  ground,  a 
built-in  approach  and  departure  system,  and  automatic 
direction  finding. 
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2.  BACKGROUND 


A  common  grid  is  3,  .system  of  readily  measurable 
coordinates  and  time  signals,  which  form  a  basic  reference 
system,  making  it  possible  for  all  users  to  make  accurate, 
real-time  ^terminations  of  present  position  and  epoch 
time.  ■  F:  .ce  all  users  of  the  system  refer  to  a  common 
reference,  transfer  of  coordinates  from  one  user  to 
another  incurs  no  loss  of  accuracy  relative  to  each  other. 
Instrumentation  to  recover  epoch  time  makes  it  feasible 
to  execute  critically  timed  ta«ks.  Lines  of  position  estab¬ 
lished  by  Loran  are. very  nearly  stationary  in  space,  in¬ 
dependent  of  weather  conditions,  diurnal  effects,  foliage 
‘  density,  terrain  roughness,  or  height  above  terrain.  This 
is  not  to  say  that  environmental  conditions  present  no  prob¬ 
lem.  Inadequate  user  equipment  or  overwhelming  noise 
may  make  it  difficult  to  obtain  a  statistically  adequate  mea¬ 
surement.  Precipitation  static  may  cause  poor  signal-to- 
noise  ratio  at  an  aircraft  antenna.  How  well  the  user  can 
make  measurements  under  difficult  environmental  condi¬ 
tions:  is  critically  dependent  on  the  design  of  system  com¬ 
ponents.  1  However,  to  the  extent  that  the  user  can  make 
the  necessary  measurement,  the  Loran  lines  of  position 
exist  as  a  stationary  grid. 

Figure  1  portrays  the  common  grid  system.  The 
Loran  service  provides  time  differences  TDA  and  TDB, 
and  the  pressure  transducer  provides  the  pressure  P. 
Clearly,  if  any  of  the  coordinates  of  two  aircraft  are 
different  simultaneously,  the  aircraft  are  not  in  the  same 
place.  Conversely,  if  one  aircraft  guides  to  a  given  set 
of  coordinates,,  such  as  TDA,  TDB, and  Pg  for  the  touch¬ 
down  point  on  a  runway,  it  will  indeed  be  there  when  its 
own  coordinates  are  the  same  as  those  of  the  touchdown 
point.  This  is  a  fundamental  principle.  All  navigation, 
horizontal  or  vertical,  can  be  done  in  a  relative  sense  in 
a  common  grid  system.  The  result  id  a  simple,  opera- 

Preceding  page  blank 
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tionally  useful  navigation  technique  of  sufficient  accuracy 
to  satisfy  the  need  in  all  phases  of  liight  with  the  possible 
exception  of  a  soft  touchdown  under  instrument  meteorologi¬ 
cal  conditions. 

Loran-C  radiates,  from  three  or  more  transmitters, 
pulsed  electrical  signals  that  are  very  accurately  timed 
and  coded  for  identification.  These  signals  are  propagated 
both  as  groundwaves  and  as  skywaves  reflected  from  the 
ionosphere.  In  the  service  area,  the  time  relationships 
that  exist  among  groups  of  signals  depend  on  where  they 
are  observed,  and  therefore  can  be  used  as  a  measurement 
of  horizontal  position. 

The  standard  measurement  technique  consists  of 
measuring  the  times  of  arrival  of  signals  from  two  trans¬ 
mitters  A  and  B  relative  to  the  time  of  arrival  of  the  signal 
from  a  third  transmitter,  M.  These  two  time  differences 
TDA  and  TDB  identify  two  hyperbolas  that  intersect  at  the 
receiving  antenna,  as  in  Fig.  2.  The  geometry  is  tied  to 
the  ground  in  the  groundwave  service  area  by  taking  into 
account  the  location  of  the  transmitters  and  the  propaga¬ 
tion  velocity.  For  skywave  service,  the  height  of  the  iono¬ 
sphere  also  must  be  taken  into  account 

If  all  the  Loran-C  transmitters  are  controlled  by 
Cesium  time  standards  and  synchronized,  it  is  expected 
that  stationary  users  of  the  signal  within  the  groundwave 
service  area  will  be  able  to  synchronize  time  to  2  ps 
(3-sigma)  and  moving  users  to  the  order  of  5^s. 

The  precision  of  position  location  available  from 
Loran  service  is  dependent  on  the  precision  of  the  time 
differences  as  they  exist  in  the  service  area.  All  Loran-C 
chains  are  monitored  at  a  site  in  the  service  area.  The 
function  of  the  monitor  is  to  compare  measured  values  of 
time  differences  to  standard  values  established  for  the 
site  and  to  use  the  differences  to  control  the  chain.  The 
best  example  of  this  process  in  operation  is  in  Southeast 
Asia  (SEA).  Analysis  of  extensive  data  recorded  24  hours 
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a  day  during  summer,  winter,  wet  and  dry  seasons  at  the 
monitor  side  at  Udorn,  Thailand,  shows  that  the  mean  of 
the  time  differences  were  controlled  within  24  nsec  of  the 
standard  control  value,  with  a  standard  deviation  less 
than  32  nsec  (see  Appendix  A).  From  the  data  of  Appendix 
A,  the  precision  of  location  on  a  geometric  dilution  con¬ 
tour  of  1  m/nsec  fell  between  8  and  40  meters  CEP.  Im¬ 
proved  monitoring  systems  planned  for  future  Loran  chains 
show  promise  of  control  substantially  better  than  is  now 
possible  by  tlie  manual  methods  in  use  in  SEA. 

Atmospheric  static  pressure  is  used  universally  as 
a  measure  of  aircraft  altitude.  The  so-called  pressure 
altitude  of  an  aircraft  is  the  ambient  static  pressure  stated 
in  terms  of  feet  derived  from  the  relationships  between 
static  pressure  and  geometric  altitude  defined  by  the  ICAO 
Standard  Atmosphere.  The  static  pressure  at  s^a  level 
for  this  atmosphere  is  29.  92  inches  of  mercury.  If  true 
geometric  altitude  is  desired,  it  must  be  estimated  by 
applying  corrections  to  the  pressure  altitude,  which  are 
themselves  estimates  of  the  differences  between  the  stan¬ 
dard  atmosphere  and  the  atmosphere  that  currently  sup¬ 
ports  the  aircraft. 
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The  problem  is  to  provide  a  simple,  reliable,  and 
inexpensive  means  for  computing  accurate  guidance  from 
Point  A  to  Point  B  along  any  path,  defined  either  in  ad¬ 
vance  (fixed  course)  or  dynamically  in  flight  (homing); 
and  to  produce  accurate  indications  of  arrival  at  specific 
points.  The  required  accuracy  of  guidance  depends  on 
the  function  to  be  performed.  It  is  commonly  stated,  for 
example,  that  the  accuracy  required  for  area  navigation 
is  far  less  critical  than  that  required  for  approaches  or 
departures.  This  is  basically  true  and  could  be  used  as 
a  basis  for  accuracy  criteria  in  any  system  in  which  area 
navigation  aids  are  different  from  terminal  area  aids. 
Suppose,  however,  that  a  single  aid  could  be  deployed, 
which  could  serve  the  critical  accuracy  requirements  of 
the  terminal  area  and  could  be  available  to  serve  the  needs 
of  area  navigation.  Such  a  system  has  great  utility  and 
economic  benefits,  since  it  reverses  the  trend  toward 
proliferation  of  navigation  aids  while  at  the  same  time 
makes  it  possible  to  execute  high-accuracy  functions  any¬ 
where.  Included  among  these  functions  are  aerial  de¬ 
livery,  aerial  survey,  traffic  control,  collision  avoidance, 
and  approach  and  departure  at  all  classes  of  airfields. 

The  next  section  describes  such  a  system.  Its  applica¬ 
bility  to  traffic  control,  collision  avoidance,  and  nonair- 
borne  uses  are  covered  in  Ref.  1. 


Preceding  page  blank 
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4.  SYSTEM  DESCRIPTION 


The  system  will  be  described  in  terms  of  the  air¬ 
craft  system  and  its  supporting  ground  system. 


AIRCRAFT  SYSTEM 

The  aircraft  system  is  shown  as  a  functional  block 
diagram  in  Fig.  3. 

Initialization 


The  data  required  for  initialization  of  the  system 
are  of  two  kinds,  shown  in  Fig.  3  as  "initial  Conditions" 
and  "Flight  Program.  "  "initial  Conditions"  consists  of 
the  numerical  values  of  a  list  of  constants  that  apply  for 
the  duration  of  a  single  Light.  These  constants  are  listed 
in  Table  1. 

The  second  kind  of  data,  "Flight  Program, "  are 
shown  in  Table  2.  This  is  organized  chronologically  in 
such  a  way  that  stepping  down  the  table  one  line  at  a  time 
provides  all  the  data  necessary  for  the  current  leg  of  the 
flight.  Looking  one  step  ahead  provides  the  data  neces¬ 
sary  to  make  the  transition  from  one  leg  to  the  next.  All 
of  the  data  of  Table  1  could  be  incorporated  in  Table  2  if 
desired,  but  this  would  result  in  substantial  repetition. 

Table  2  is  organized  to  permit  the  flight  to  pass 
through  any  number  of  Loran  triads  and  chains.  It  also 
allows  any  three  transmitters  to  be  designated  a  triad. 
Most  Loran  receivers  operate  only  on  one  group  repeti¬ 
tion  period  at  a  time.  Therefore,  with  these  receivers, 
all  three  transmitters  would  have  to  be  members  of  one 
chain.  Also,  most  receivers  require  the  master  trans¬ 
mitter  to  be  a  member  of  e?ch  triad.  Neither  of  these  re¬ 
strictions  is  necessary  in  receivers  of  the  future. 

Preceding  page  blank 
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TABLE  1 

CONSTANTS  FOR  INITIAL  CONDITIONS 

a.  Loran  Receiver 

P,  The  group  repetition  period  of  the  I**1  chain; 

D}  The  coding  delays  for  the  transmitter  on  the 

Ith  chain,  J  =  1.  2,  3  ...  M. 

b.  Conversion  to  Loran  Rectangular  Coordinates 

j,  Yj.  j  The  Loran  rectangular  coordinates  of  the  Jth 
1  ’  transmitter  in  the  Ith  chain,  1  =  1,  2  ...  N; 

J  =  0,  1,  2  . . .  M,  where  J  =  0  designates  the 
master; 

V  The  propagation  velocity  of  the  Loran  ground 

wave. 

c.  Balnstic  Data 

There  are  two  alternative  initial  conditions  depending  on  the 
choice  of  ballistics  algorithm  embodied  in  the  avionic  com¬ 
puter. 

c-1.  The  --simplest  algorithm  designed  for  horizontal  releases  re¬ 
quires  the  following  data; 

TRq  ,  TFq  The  TRail  and  Time  of  Fall  of  the  store  when 
released  under  the  planned  conditions  of  re¬ 
lease; 

TRV,  TFV  The  partial  derivatives  of  TR  and  TF  with  re¬ 
spect  to  Velocity  pressure  taken  at  the  planned 
true  air  speed  at  release; 

TRH,  TFH  The  partial  derivatives  of  TR  and  TF  with  re¬ 
spect  to  Height  pressure  at  the  planned  altitude 
of  release; 

VPq  The  desired  Velocity  Pressure  at  release 

^  The  compass  variation  plus  deviation. 
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TABLE  1  (cont'd) 


c-2.  The  more  complex  algorithm  requires  the  following  data: 


d. 


e. 


bik 

$ 

Approach 

tan»4 


General 

a 

TCR 

max 


tan  o 


max 


The  ballistic  index  number  that  identifies 
store  K; 

The  compass  variation  plus  deviation. 


The  tangent  of  the  glide  slope  for  approach; 

The  slope  of  pressure  versus  altitude  at  sea 
level; 

The  correction  factor  to  obtain  the  slope  of 
pressure  versus  altitude  at  altitudes  above 
sea  level. 


Nominal  speed  of  sound; 

The  maximum  horizontal  angular  rate  of  turn 
for  the  aircraft; 

The  tangent  of  the  maximum  bank  angle  of 
the  aircraft. 
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Table  2  also  shows  an  example  of  the  listings  re¬ 
quired  in  three  columns  for  a  particular  flight.  In  the 
"Mode"  column,  the  flight  starts  in  the  area  navigation 
mode  and  proceeds  in  that  mode  until  waypoint  j.  Here  is 
entered  the  designator  2  for  delivery  mode,  indicating  that 
the  leg  of  the  flight  approaching  waypoint  3  is  a  delivery 
leg,  and  guidance  is  to  be  provided  to  an  air  release  point 
such  that  delivery  will  be  made  to  waypoint  3'  associated 
with  target  pressure  GP^.  In  the  "BI"  column,  a  ballistic 
index  number  is  listed  also  against  waypoint  3,  indicating 
the  specific  store  to  be  released,  assuming  the  more  com¬ 
plex  ballistic  algorithm  is  used  (see  Sect,  c-2  of  Table  1). 
If  the  simpler  algorithm  (Sect,  c-1  of  Table  1)  is  used,  the 
"BI"  column  is  not  needed.  If  the  designator  3  for  recon¬ 
naissance  were  entered,  guidance  would  be  to  waypoint  3 
and,  on  arrival,  a  signal  for  gathering  data  would  be  pro¬ 
duced.  In  the  "Location  (Next  Triad)"  column,  a  pair  of 
coordinates  are  shown  against  waypoint  1 ,  which  are  the 
coordinates  in  the  next  triad,  which  corresponds  to  the 
same  point  as  the  coordinates  shown  in  the  "Location 
(Present  Triad)"  column.  When  this  double  listing  of 
coordinates  occurs,  it  indicates  that  the  next  leg  of  the 
flight  beyond  waypoint  1  is  to  be  flown  on  a  course  from 
X[,  Yj  to  X2J  Y2  in  the  new  triad.  Finally,  in  the  "Mode" 
column,  designator  4  is  listed  for  approach.  Asa  conse¬ 
quence,  the  last  leg  is  flown  on  a  glide  path  to  waypoint 
n  and  ground  pressure  GPn. 

Loran  Receiver 


Refer  now,  in  Fig.  3,  to  the  block  "Loran  Receiver.  " 
The  initial  conditions  for  the  receiver,  given  in  Table  1,  are 
set  to  the  values  required  by  the  flight  program.  Table  2, 
and  the  receiver  produces  the  Loran  coordinates,  TDA  and 
TDB,  of  present  position  and  the  time  rates  of  change  of 
TDA  and  TDB.  (Details  of  how  Loran  receivers  produce 
hyperbolic  coordinates  and  their  time  rates  are  beyond  the 
scope  of  the  pre  sent  effort. ) 
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Loran  Rec  angular  Coordinates 

The  next  block  is  "Loran  Rectangular  Coordinates.  " 
The  function  of  this  block  is  to  convert  the  Loran  hyper¬ 
bolic  coordinates  to  Loran  rectangular  coordinates.  This 
is  done  as  follows: 


X__,  Y.. 
M  M 

XA.  Y. 
A  A 

XB'  YB 

XR*  YR 


TDA 

TDB 

DA 


be  the  coordinates  of  the  Master, 

be  the  coordinates  of  Slave  A, 

be  the  coordinates  of  Slave  B, 

be  the  coordinates  of  the  receiver  corre¬ 
sponding  to  TDA,  TDB, 

be  the  time  difference  A, 

be  the  time  difference  B, 

be  the  absolute  emission  delay  for 
Slave  A,  and 

be  the  absolute  emission  delay  for 
Slave  B. 


Then,  for  the  current  chain  and  transmitters  identified  by 
Table  2, 


x„. 

Y__ 

=  X  T 

M 

M 

I,J 

XA- 

ya 

■ XI.  J 

XB' 

yb 

-  XI,  J 

DA 

=  DI,J 

DB 

=  DI,J 

and  the  following  expressions  (us  ing  self-explanatory 
dummy  variables  for  clarity)  produce  the  values  of 
xr,  Yj^  corresponding  to  TDA  and  TDB  (see  Appendix  B): 
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MAMTUNO 


A 

=  Y  -  Y 

M  XA 

B 

ii 

& 

C 

=  XM  "  XA 

D 

=  XM"XB 

A 

=  A-D  -  B-C 

=  V* (TDA  -  DA) 

=  V*  (TDB  -  DB) 

E 

=  d*ma 

F 

=  c*mb 

Of 

=  0.  5(|4a*  E  -  F)/A 

=  (E  -  F)/A 

G 

=  a-mb 

H 

=  b*ua 

'/ 

=  0.  5(mb'G  -ua*H)/A 

6 

=  (G  -  H)/A 

$ 

=  (YM-«)2  +  (XM-y)2 

z 

V 

H 

1 

1 

cO 

+ 

1 

aa. 

it 

=  <  We2  -  4£t? 

2n 

Negative  solutions  for  0^  are  not  real.  Therefore,  these 
can  be  eliminated  by  a  simple  sign  check.  However,  there 
are  regions  in  the  service  area  where  two  positive  values 
are  obtained  for  0j^.  In  this  case,  a  priori  knowledge 
must  be  used  to  select  the  correct  solution  as  follows: 

Let 


XQ ,  Yq  be  the  coordinates  of  the  receiver  position 
at  the  beginning  of  a  flight  (see  Table  2>. 
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Let 


®M " (8M*2 ' 


if  previous  0  *s  had  been  (0.„), . 

M  Ml 


Then  let 


WM*0  ~  0M  ’ 


and  solve  for 

XR  ">'+6‘eM- 


yr=°+^9m- 


•  The  velocity  components  of  the  receiver,  Xj>  an<* 

Yr  arc  obtained  as  follows: 


Let 


a  = 


XR  -  XA  XR-XM 


b  = 


c  = 


d  = 


0--  +  Ma 

M  A 

3  • 

M 

yr  "  ya 

yr  "  ym 

8m  +  ma 

0M 

XR-XB 

X_  - 
R  M 

®M  +  UB 

0,, 

M 

yr‘yb 

yr  ‘  ym 

0m~  ub 

®M  ’ 
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and 


AA  = 


ad  -  be 


BB  = 


-b 


ad  -  be  * 


CC  = 


-a 


ad  -  be  * 


DD  = 


ab  -  be  * 


Then 


X„  =  V(AA*  TDA  +  BB*  TDB) 

XV 

Y0  =  V(CC*  TDA  +  DD*  TDB) 

XV 


where 

•  • 

TDA  and  TDB  are  the  time  rates  of  change  of  TDA 
and  TDB. 


The  logic  given  above  for  selecting  the  sign  of  the 
square  root  in  '  ^  solution  for  0^j  is  involved  with  an  in¬ 
determinacy  whenever  =  (6m^2*  This  equality  only 

exists  along  base  line  extensions.  (In  a  triad,  there  are 
three  base  lines  and  six  baseline  extensions.  An  extension 
is  a  line  drawn  from  a  transmitter  in  a  direction  opposite 
to  the  direction  of  any  other  transmitter. )  If  the  logic 
given  above  is  used  when  (6^)1  anc*  are  e<lual»  the 

indeterminacy  is  resolved  for  cases  where  the  baseline 
extension  is  being  crossed  at  a  reasonable  angle.  How¬ 
ever,  if  a  sufficient  length  of  time  is  spent  continuously 
in  close  proximity  to  the  baseline  extension,  the  naviga¬ 
tion  computer  will  lose  track  of  which  way  it  should  leave 
the  baseline  extension.  This  problem  can  be  resolved 
only  with  additional  a  priori  ii  formation.  The  logic  re¬ 
quired  to  resolve  this  complication  will  not  be  developed 
in  this  paper  because  Loran  along  baseline  extensions  is 
known  to  be  poor  for  other  reasons.  The  preferred 
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resolution  of  the  problem  is  to  place  operational  constraints 
on  aircraft  and  Loran  chain  design  so  that  aircraft  never 
have  to  fly  parallel  and  close  to  baseline  extensions. 

Navigation 

Refer  now  to  the  block  labeled  "Navigation”  in  Fig. 

3.  The  inputs  from  "Loran  Rectangular  Coordinates"  are 
a  continuous  succession  of  the  receiver  coordinates  Xr 
and  Yr  and  its  velocity  components  X^  and  Yr.  Another 
input,  prior  to  the  start  of  each  leg  independent  of  mode, 
is  the  waypoint  information  from  the  "Flight  Program"  of 
Table  2.  In  the  case  of  the  delivery  mode,  the  upcoming 
waypoint  is  adjusted  to  the  coordinates  of  the  release  point. 


a.  Fixed  Course.  The  navigation  problem  is 
"getting  from  waypoint  j  to  waypoint  j  -t- 1  *’  for  any  value 
of  j  from  0  to  n  -  1.  Figure  4  displays  *his  general  situa¬ 
tion  with  the  desired  guidance  data.  Cross  Track  Error 
(CTE),  and  Mong  Track  Distance  to  go  (ATD).  This  navi¬ 
gation  mode  is  called  "Fixed  Course. " 


Let  us  denote  Xj+j  -  Xj 
AYj,  as  in  Fig.  4.  The  true  cc 
ana  the  following  relations  hold 


-  Xj  as  AXj  and  Yj+j  -  Yj 
e  true  course  from  j  to  i+1  is 


TCj, 


sinTC. 

3 


I  2  2 
\/£X.  +  AY 
'  3  3 


cosTC.  = 
3 


2  2 
AX  +  AY 
3  3 


Now,  the  current  position  of  the  receiver  is  Xjj,  Y^.so 
that 


Ae  =  X  -  X_  ; 
3+1  3+1  R 
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NORTH 


WAYPOINT 

i  +  1 


Fig.  4  HORIZONTAL  FIXED  COURSE  GEOMETRY 
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An  =  Y  -  Y  : 
•  j+1  j+1  R  ' 


where  Ae  and  An.+  1  are  the  east  a^d  north  positions  of 
the  receiver  with  respect  to  the  j  +  lwaypoint.  The  steer¬ 
ing  signals  are  now  given  by 


CTE . , ,  =  Ae . , ,  cosTC .  -  An. ,  .  sinTC .  , 
J+1  J+1  3  3+1  3 


ATD.,  ,  =  Ae.,,  sinTC.  +  An.it  cosTC.  , 
3+1  3+1  3  3+1  3 


CTER., ,  =  -XL,  cosTC.  +  Y„  sinTC. 
3+1  Tt  3  R  3 


ATDR  =  -XD  sinTC.  -  Y-,  cosTC. 
J+1  R  3  R  3 


where  CTER  and  ATDR  are  time  rates  of  change  of  CTE 
and  ATD.  If  required. 


HE  =  tan 


-1  CTE 
ATD  1 


where  HE  is  the  heading  error. 


f.  i 


b.  Homing.  Fixed  course  guidance  has  been  de¬ 
scribed  above.  Another  option  is  "Homing"  guidance.  In 
this  case.  Fig.  5  pertains,  and 


sinTC. 

J 


cosTC . 
J 


/  2  2 

/Ae*  +  An.^, 
3+1  J+1 


/  2  2 

Ki + AVi 
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The  ADF  function  (Automatic  Direction  Finding)  is  fulfilled 
by  taking  ™ 


TC.  =  sin 
J 


-1 


Ae 


3+1 


4 


A Vi + 


The  steering  signals  are  now  given  by: 


CTBj+l  =  0  , 

ATI V =  AVi 

CTERj+l  *  -*R 

ATDV  = 


sinTC.  +  An.  ,  cosTC.  , 
3  3+1  3 

• 

cosTC .  +  sinTC .  , 

3  K  ] 

sinTC.  -  Y„  cosTC.  . 

3  R  3 


c.  Turns.  If  a  transition  from  one  flight  leg  to 
the  next  is  not  initiated  until  the  upcoming  waypoint  is 
reached,  the  turn  will  result  in  considerable  overshoot. 
The  point  to  initiate  the  turn  in  advance  of  the  upcoming 
waypoint  is  determined  as  follows:  Aircraft  are  often 
steered  by  executing  standard  rate  turns  unless  the  de¬ 
mands  on  the  aircraft  require  bank  angles  greater  than  a 
stated  maximum,  in  which  case  the  maximum  bank  angle 
is  executed.  Thus, 


TCR  *  TCR 

max 

and 


where  TCR  is  the  time  rate  of  change  of  course  and  <p  is 
the  angle  of  bank  in  a  steady-state  turn.  First  find 

TAS  =  Ma , 

_  TAS 
1  "  TCR 

max 
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Then  find 


r  -  TAS 
2  g  tan  <p 


max 


where 

TAS  is  the  True  Air  Speed, 
a  is  the  speed  of  sound, 

R  is  the  radius  of  turn,  and 
g  is  acceleration  of  gravity. 

If  Rj  £  R2>  execute  the  turn  at  a  point  determined  by  Rj. 
If  R^  >  R2,  execute  the  turn  at  a  point  determined  by  R2. 
Now  let  R  be  either  Rj  or  R2,  depending  on  the  above 
conditions,  and  refer  to  Fig.  6. 

Let 


and 


ATC  =  TC.,, 
3+1 


9 


AATD  =  R  tan 


ATC 
2  * 


Depending  on  the  sign  of  ATC,  a  right  or  left  turn  is  exe¬ 
cuted  at  the  Turning  Point,  TP,  that  is,  when  ATD  = 
AATD.  There  are  four  options  for  making  the  turn.  The 
first  consists  of  flying  the  standard  turn  associated  with 
either  TCjj^x  or  <pm  until  such  time  as  the  next  leg  of 
flight  progresses  to  AATD  from  the  starting  waypoint,  at 
which  time  the  standard  turn  is  terminated  and  guidance 
resumed  according  to  CTE  indications.  The  second  op¬ 
tion  is  to  change  at  TP  from  Fixed  Course  to  Homing  on 
waypoint  j+2,  then  change  back  to  Fixed  Course  when 


(CTER 


\ 

j+2'Homing 


=  0  ±  (CTER) 


Tolerance 
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Having  passed  waypoint  j+1,  the  Fixed  Course  guidance 
will  be  along  the  leg  j+1  to  j+2.  The  thirc  option  is  to 
change  to  Fixed  Course  guidance  along  leg  j+1  to  j+2  at  TP, 
at  which  time  substantial  CTE  and  CTER  will  exist,  indi¬ 
cating  the  correct  guidance  to  the  desired  course.  This 
technique  can  produce  very  accurate  ground  tracks  in  a 
turn  by  spacing  the  waypoints  very  closely  during  the  turn, 
and  it  is  the  simplest  to  execute.  If  guidance  is  required 
along  a  continuous  curved  path,  a  fourth  option  would  be  to 
compute  on-board  the  geometry  of  the  curve  desired  for 
the  turn. 

d.  Vertical  Guidance.  Another  input  to  the  "Navi¬ 
gation"  block  is  atmospheric  static  pressure,  i.  e. ,  the 
height  pressure  of  the  aircraft  HP,  and  the  height  pres¬ 
sure  time  rate  HPR.  Also,  an  input  from  the  Flight  Pro¬ 
gram  is  the  height  pressure  associated  with  tue  waypoints. 
There  are  basically  two  ways  to  handle  vertical  guidance 
in  the  area  navigation  mode.  The  first  way  is  to  produce 
continuous  altitude  indications  analogous  to  the  continuous 
horizontal  guidance  described  above.  The  second  way  is 
to  produce  incremental  vertical  indications  between  way- 
points. 

Figure  7  displays  the  general  situation  with  the 
desired  vertical  guidance  data.  Height  Pressure  Error 
(HPE)  and  Height  Pressure  Error  Rate  (HPER).  Figure  7 
is  the  situation  in  a  vertical  pla .  e  pa  ssed  through  the  way- 
points  of  Fig.  4.  Let  us  denote  HP-j+j  -  HPj  as  AHPj. 

The  desired  climb  path  angle  from  j  to  j+1  is  y,  and  Uie 
following  relation  holds: 

AHP. 

K  tany  = 


Now,  the  current  position  of  the  pitot  is  HP,  so  that 

HPE  =  HP  -  ATD  K  tany  -  HP  . 

3+1  3+1  3+l 
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Fig.  7  VERTICAL  GUIDANCE  GEOMETRY 


TMC  JOHN*  HOPKINS  UNtVtfttlTY 

APPLIED  PHYSICS  LABORATORY 
suvt*  •»— —  mavtuno 


Similarly 

HPER .  ,  =  ATDR  K  tany  -  HPR. 

3+1  3+1  r 

This  relation  holds  for  all  j,  and  does  not  require  any 
trigonometric  routines.  Another  set  of  outputs  of  the 
"Navigation"  block  is  height  pressure  error  and  height 
pressure  error  rate.  In  flight  planning,  restrictions 
must  be  placed  on  y  so  that  the  demands  on  the  aircraft 
do  not  exceed  its  capability. 

Referring  again  to  Fig.  7,  the  incremental  verti 
cal  guidance  techniques  consist  of  letting 


HPE  =  HP_,  -  HP  . 

3+1  3+1 

Then,  depending  on  the  sign  of  HPE,  a  standard  climb  or 
glide  is  executed  until  HPEj+j  =  0.  At  this  point,  the 
climb  or  glide  is  terminated  and  level  flight  achieved, 
such  that  HPE.+1  continues  to  be  zero.  In  this  case,  dur¬ 
ing  flight  planning  the  waypoints  could  be  placed  so  close 
together  that  the  aircraft  cannot  reach  the  desired  altitude 
by  the  time  it  begins  the  next  leg.  This  is  of  no  conse¬ 
quence  unless  being  at  altitude  throughout  the  leg  j+1  to 
j+2  is  important  to  the  mission. 

e.  Arrival  Time.  Under  some  circumstances  it 
may  be  important  to  arrive  at  waypoints  at  specified  times. 
If  this  is  the  case,  the  required  ground  speed  along  each 
leg  of  the  flight  is  determined  as  follows: 


ATDRV  *  T 


ATD 

— -  ATDV  • 


where 


T 

3+1 


R 


3+1 


R 


is  the  required  Time  Of  Arrival  (TOA)  at 
waypoint  j+1, 

is  time  at  the  receiver,  and 


ATDRE  is  the  along-track  distance  rate  error. 
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During  mission  planning,  the  time  of  arrival  at  each  way- 
point  must  be  related  to  the  length  of  the  leg  so  that,  the 
speed  capability  of  the  aircraft  is  not  exceeded,  considering 
the  influence  of  winds  aloft. 

Air  Data 

The  vertical  positioning  of  the  aircraft  and  the  hori¬ 
zontal  location  of  an  air  release  point  are  dependent  on  the 
accuracy  of  the  air  data.  Two  pressure  measurements  are 
required,  ram  and  static.  The  measurement  of  true  ram 
pressure  is  not  difficult  since  a  true  value  is  obtained  when 
the  air  is  suddenly  brought  to  rest  anywhere  outside  the 
boundary  layer.  However,  very  accurate  measurements  of 
true  static  pressure  are  virtually  impossible  because  of  the 
potential  flow  interference  effects.  These  interference  ef¬ 
fects  must  be  accounted  for  through  accurate  calibration  of 
the  static  pressure  measurements  for  the  specific  location 
of  the  static  orifices,  for  the  operational  speed  range,  and 
for  operationally  important  flight  configurations,  e.  g. ,  ! 
flaps  versus  no  flaps. 

Let 

P  be  the  ram  pressure, 

A 

P  be  the  static  pressure  measurement, 

b 

AP  be  the  static  pressure  interference  effect, 

HP  be  the  true  static  pressure  or  Height 
Pressure, 

HPR  be  the  Height  JPressure  time  Rate  of  change 

"  i 

VP  be  the  velocity  pressure, 

t  be  the  time, 

M  be  the  Mach  number,  and 

y  be  the  ratio  of  the  specific  heats  for  air. 

Then,  the  desired  outputs  of  the  "Air  Data"  block  are 
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HPR 

VP 

I 


There  are  two  options.  If  the  geometric  altitude  of  a 
•  specific  type  .of  aircraft  does  not  have  to  be  accurately 
controlled,  the  output  of  the  "Air  Data"  block  can  be  "raw" 
data,  i.  e. ,  AP  =  0.  In  this  case,  the  effects  of  AP  can  be 
accounted  for  in  the  ground  system.  Under  the  other  op¬ 
tion,  the  "Air  Data"  block  output  is  as  accurate  as  the  caii- 
■  bration,  i.  e. ,  AP  is  calculated  as  a  function  of  Mach  num¬ 
ber  dnd  aircraft  trim  condition.  In  this  case,  thp  ground 
system  ignores  the  effects  of  AP. 

:  :  i 

Air  Release  Point 

The  path  of  an  unguided  object  released  from  an  air¬ 
craft  is  dependent  on  its  aerodynamic  characteristics  and 
the  initial  conditions  at  release.  To  ladnch  a  specific  ob¬ 
ject  on  a  predicted  path,  the  initial  conditions  must  be  con- 
,  trolled  to  specific  values;  and  to  simplify  this  task,  releases 
are  usually  made  tinder  steady-state  conditions.  The  fol¬ 
lowing  computation  of  air  release  point  applies  to  steady 
conditions  of  true  air  speed  and  vertical  velocity,  i.  e. ,  the 
only  acceleration  acting  on  the  aircraft  is  gravity. 

Let 

i 

TAS  be  the  T^rue  _Air  Speed, 

TH  be  the  True  Heading,; 

TP  be  the  Time  of  Fall, 

GS  be  the  Ground  Speed, 

*  TR  ,  be  the  TRail, 

BA  be  the  Ballistic  Altitude, 

DR  be  the  DRift  Angle, 

i 
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TC  be  the  True  Course, 

CT  be  the  Cross  Trail, 

AR  be  the  Actual  Range, 

W  be  the  Vertical  Velocity, 

BI  be  the  Ballistic  Index  number, 

BR  be  the  Ballistic  Range,  and 

if)  be  the  compass  variation  plus  deviation. 

Then,  referring  to  Fig.  8, 

TH  =  H  +$, 

DR  =  TC  -  TH, 

CT  =  TR  sin  DR, 

AR  =  TF  •  GS  -  TR  cos  DR, 

where  H  is  the  compass  Heading.  Table  2  indicates  that 
the  target  horizontal  coordinates  are  given  as  one  of  the 
waypoints  identified  by  "2"  for  delivery  mode.  When  this 
designation  occurs,  the  guidance  for  that  leg  of  the  flight 
is  computed  by: 

CTE  =  CTE.^,  +  CT, 

K  j+l 

ATDd  =  ATD. ,,  -  AR, 

K  J+l 

where  CTEj^  and  ATDr  are  with  respect  to  the  air  Release 
point. 


As  indicated  in  Table  1,  there  are  two  options  for 
the  ballistic  solution.  The  first  requires  the  aircraft 
vertical  velocity  to  be  zero,  and  it  is  linearized  over 
limited  deviations  of  altitude  and  air  speed  from  the  values 
specified  by  Tables  1  and  2.  The  vertical  and  horizontal 
components  of  ejection  velocity  are  accounted  for  by  TRq 
and  TFq  .  The  following  equations  apply: 
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Fig.  8  BALLISTIC  SOLUTION 


A  VP  =  VPQ  -  VP, 

AHP  =  HP.l1  -  HP, 

J+l 

TR  =  TRq  +  TRV  •  A  VP  +  TRH  •  AHP, 

TF  =  TFq  +  TFV  •  A  VP  +  TFH  •  AHP, 

where  TRV,  TFV  are  the  derivatives  of  TR  and  TF  with 

respect  to  Velocity  pressure  taken  at 
the  planned  true  air  speed  at  release; 
and 

TRH,  TFH  are  the  derivatives  of  TR  and  TF  with 
respect  to  Height  pressure  taken  at 
the  planned  altitude  of  release. 

The  second  option  requires: 

BP  =  GP  -  HP, 

where  BP  is  the  Ballistic  Pressure  height  above  ground. 
Ballistic  range  and  time  of  fall  are  determined  as 

TR  =  f,  (VP,  HPR,  BP,  BI), 

1 

TF  =  f2(VP,  HPR,  BP,  BI), 

where  the  functions  and  f2  are  arbitrary  algebraic  ex¬ 
pressions  that  represent  the  ballistic  tables  for  any  object 
released  with  specified  ejection  velocity  within  the  flight 
regime  of  the  aircraft  under  steady-state  conditions.  The 
vertical  velocity  effect  can  be  accounted  for  either  as  a 
rate  of  change  of  Dressure  HPR  from  Air  Data  or  as  actual 
vertical  velocity  converted  to  HPR  in  aircraft  equipped 
with  suitable  inertial  components.  If  the  "Air  Data"  block 
of  Fig.  3  removes  the  effect  of  AP,  the  coefficients  of  the 
equations  for  TR  and  TF  reflect  "true"  values.  If  AP 
effect  is  not  removed  in  the  aircraft,  the  coefficients  must 
reflect  the  calibration,  and  then  the  equations  for  TR  and 
TF  apply  only  to  that  aircraft. 
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Glide  Slope 

Refer  to  Fig.  9  for  the  approach  geometry  where  7^ 
is  the  desired  glide  slope  angle  and  AH  is  the  desired  in¬ 
cremental  height  above  waypoint  n  located  at  the  touch¬ 
down  point.  Now  proceed  as  follows: 

AH  =  ATD  tan7.  , 

n  n  A 

HP  =  GP  +  [Sn  +  SJP^  -  GP  )]  AH 

n  n  L  0  I  S  n  n 

and 

HPE  =  HP  -  HP  , 
n  n 

where 

Sq  is  the  Slope  defined  as  the  derivative  of  atmo¬ 
spheric  pressure  with  respect  to  height  at  sea 
level, 

Sj  is  a  correction  factor  for  determining  the  Slope 
at  a  height  equivalent  to  GP,  and 

Pg  is  the  sea  level  atmospheric  pressure  of  the 
standard  atmosphere. 

If  values  of  SQ  and  Sj  are  determined  from  the  ICAO  stan¬ 
dard  atmosphere  using  inches  of  mercury  and  feet  of  alti¬ 
tude,  SQ  =  -  1. 05  x  10"  3  and  Sx  =  2.  58  x  10-5. 

Aircraft  Control 


The  output  of  the  "Navigation”  block  of  Fig.  3  is 
•Cross  Track  Error  (CTE),  Cross  Track  jError  Rate  (CTER), 
Height  ^Pressure  jError  (HPE),  and  Height  Pressure  jError 
jftate  (HPER).  These  outputs  must  be  transformed  into  con¬ 
trol  commands  to  either  an  autopilot  or  pilot.  In  the  case 
of  autopilot  control,  the  commands  are  usually  accelerations 
that  are  obtained  from  equations  of  the  form 


Cx  =  A^CTE)  +  Bj(CTER), 
c2  =  A2(HPE)  +  b2(hper>. 


where 


1,2 

Al,2andB1.2 


are  the  acceleration  commands,  and 

are  coefficients  that  are  selected  to 
produce  the  desired  performance. 
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In  the  case  of  control  by  the  pilot,  the  commands  are  dis¬ 
played  on  the  bank  and  pitch  steering  bars  of  the  Attitude 
Director  Indicator  (ADI),  and  CTE  is  displayed  by  the 
Horizontal  Situation  indicator  (HSI).  In  either  case,  the 
desired  "Aircraft  Dynamics"  are  obtained  through  the  con¬ 
trol  actuators.  The  feedback  loops  to  Antenna  location. 
Heading,  Air  Data,  and  possibly  Autopilot  characteristics 
are  closed  through  the  aerodynamics. 

In  case  of  failure  of  the  primary  navigation  system, 
back-up  navigation  aids  may  be  available.  If  they  are, 
outputs  must  be  transformed  to  the  same  type  of  flight 
control  commands  that  the  autopilot  or  cockpit  displays 
receive  from  the  primary  system.  Their  commands  would 
be  introduced  at  the  point  indicated  by  the  dashed  line  in 
Fig.  3.  The  flight  program,  initialization,  and  algebraic 
description  of  backup  systems  are  beyond  the  scope  of 
this  paper. 

Mode  Characteristic s 


When  the  aircraft  is  operating  in  any  of  the  modes, 
the  system  continuously  produces  guidance  signals  to  fly 
the  flight  profile  defined  by  the  Flight  Program,  and  signals 
are  produced  at  each  turning  point  to  sequence  the  "Flight 
Program"  to  the  next  leg  of  flight.  In  addition,  when  the 
aircraft  is  operating  in  the  delivery  mode,  a  release  signal 
is  produced,  and  pertinent  data  are  frozen  in  the  navigation 
computer  at  the  computed  air  release  point.  The  frozen 
data  are  then  recorded.  When  operating  in  the  reconnais¬ 
sance  mode,  signals  are  produced  to  initiate  and  terminate 
the  acquisition  of  data  at  a  waypoint  identified  by  the  flight 
program.  Also,  appropriate  flight  data  are  frozen  in  the 
navigation  computer  and  recorded. 


GROUND  SYSTEM 


The  ground  system  required  to  produce  flight  pro¬ 
grams  is  shown  as  a  functional  block  diagram  in  Fig.  10. 
Basic  data  required  for  target  and  waypoint  coordinates 
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Fig.  10  GROUND  SYSTEM 
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are  obtained  from  the  two  functional  blocks  "Target  System" 
and  "Cartographic  Support.  "  The  distinction  is  that  the 
"Target  System"  function  is  to  produce  the  coordinates  of 
specific  points  with  the  highest  accuracy  that  is  practicable 
to  the  aircraft  system  operating  under  critical  guidance 
conditions.  The  "Cartographic  Support"  function  is  to  pro¬ 
duce  cartographic  materials  that  allow  operators  to  deter¬ 
mine  coordinates  of  less  accuracy  than  the  "Targeting  Sys¬ 
tem"  but  of  sufficient  accuracy  for  noncritical  guidance 
conditions.  "Meteorological  Support"  provides  the  infor¬ 
mation  required  to  convert  altitude  to  height  pressure. 
"Technical  Manuals"  and  "Planning  Materials"  provide 
instructions,  ballistic  data,  and  tabulation  forms  for  spe¬ 
cifying  the  flight  profile  and  other  flight  data. 

Target  System 

The  basic  aircraft  reference  system  depends  on 
Loran  time  differences  and  height  pressure.  Basic  tar¬ 
geting  systems  depend  on  geodetic  coordinates  and  geo¬ 
metric  altitude  above  sea  level.  Therefore,  a  coordinate 
transformation  is  required  to  properly  interface  the  two 
systems. 

In  the  past  the  classical  transformation  has  been 
used,  which  depends  for  accuracy  on  accurate  determina¬ 
tions  of  geodetic  distances  and  propagation  velocities. 
Geodetic  distances  depend  in  turn  on  mathematical  models 
of  the  earth,  and  propagation  velocities  depend  on  de¬ 
tailed  knowledge  of  geology  and  its  influence  on  earth  con¬ 
ductivity.  Experience  has  shown  that,  in  spite  of  its  com¬ 
plexity,  the  accuracy  of  the  classical  transformation  does 
not  fully  exploit  the  accuracy  potential  of  Loran-C  without 
empirical  correction  factors  derived  from  calibration 
data.  Analysis  of  errors  made  by  the  classical  transfor¬ 
mation  shows  that  the  significant  errors  are  attributable  to 
imperfections  of  the  propagation  model.  In  view  of  these 
imperfections  and  the  necessity  to  calibrate,  it  is  not 
necessary  to  continue  the  use  of  sophisticated  geodetic  and 
propagation  models  for  purposes  other  than  research.  The 
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desired  transformation  from  hyperbolic  to  geodetic  coor¬ 
dinates,  or  the  inverse,  is  fully  realized  by  the  simplest 
mathematical  relationship  that  preserves  the  accuracy  in¬ 
herent  in  the  calibration.  The  following  expressions  capi¬ 
talize  on  this  concept  and  greatly  reduce  the  computer 
burden  associated  with  the  transformation  function. 

Assume  that  calibration  data  are  gathered  at  dis¬ 
tinctive  geographic  features  whose  geodetic  coordinates 
are  accurately  known.  These  can  be  any  standard  set, 
the  most  common  being  Longitude/ Latitude,  and  "UTM" 
(Universal  Transverse  Mercator).  Using  UTM  as  an  ex¬ 
ample,  the  calibration  data  are  a  collection  of  two  pairs 
of  coordinates,  UTM  and  hyperbolic,  both  of  which  apply  to 
the  same  geographic  feature.  The  UTM  coordinates  E 
and  N  are  known  from  prior  mapping  and  the  hyperbolic 
coordinates  TDA  and  TDB  are  measured.  Now  let 

TDA  =  An  +  A,  E  +  A0  N  +  A.  E-  N  +  A.E2  +  A.N2  , 
0  1  2  3  4  5 

TDB  =  B  +  B.  E  +  B„  N  +  B„  E*  N  +  B  .E2  +  BCN2  , 

U  1  Z  o  4  O 

and  find  the  values  of  the  A  and  B  coefficients  by  a  least 
square  curve  fit  process  to  the  calibration  data.  Substitu¬ 
tion  of  target  coordinates  E.p  and  N-p  in  these  equations 
produces  values  of  TDA  and  TDB  with  an  accuracy  that  is 
at  least  as  good  as  the  calibration  data.  In  cases  where 
redundant  data  are  obtained,  the  accuracy  of  TDA  and  TDB 
enjoys  a  statistical  improvement. 

The  targeting  system  also  produces  target  altitude 
as  geometric  height  above  sea  level.  Conversion  to  height 
pressure  is  not  made  by  the  "Targeting  System"  because 
this  conversion  is  most  accurate  if  made  with  respect  to 
the  atmospheric  pressure  distribution  prevailing  during 
the  time  of  the  flight  in  the  target  area.  Due  to  temporal 
changes  in  the  atmosphere,  this  is  best  deferred  as  long 
as  possible  prior  to  takeoff  or,  if  the  system  is  so  de¬ 
signed  in  advance,  the  conversion  could  be  made  after 
takeoff  and  communicated  to  the  aircraft. 
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Cartographic  Support 

For  general  operational  utility,  the  inverse  of  the 
above  equations  for  TDA  and  TDB  can  be  used  to  draw  a 
grid  of  TDA  and  TDB  on  the  standard  maps  normally  used 
for  laying  out  flight  plans.  The  accuracy  of  coordinate 
transformation  accomplished  in  this  manner  is  degraded 
somewhat  by  the  necessity  to  visually  recall  the  data  from 
the  maps,  and  by  the  fact  that  maps  never  represent  the 
current  geographic  situation.  However,  such  graphics 
satisfactorily  fulfill  the  need  for  noncritical  v.aypoint  loca¬ 
tions  and  general  flight  planning. 

Meteorological  Support 

The  function  of  "Meteorological  Support"  is  to  sup¬ 
ply  the  atmospheric  data  required  to  convert  altitude  to 
height  pressure.  Consistent  with  the  equations  shown  be¬ 
low  under  "Height  Pressure, "  the  data  required  are  the 
sea  level  pressure  in  the  target  areas  (or  any  area  w  ere 
vertical  aircraft  control  is  critical)  and  an  indication  of 
the  prevailing  equations  of  state  for  the  area  of  operations; 
for  example,  the  equations  for  any  one  of  the  standard  at¬ 
mospheres  of  Ref.  2. 

Technical  Manuals 


The  required  manuals  are  the  standard  types  that 
describe  the  system  and  how  to  operate  it.  The  following 
paragraphs  describe  that  part  of  the  manuals'  content  that 
is  not  obvious  from  other  sections  of  this  report. 

In  the  case  of  the  simpler  ballistic  solution  (see 
Table  1),  the  manuals  must  include  ballistic  tables  that 
list  Trail  and  Time  of  Fall  versus  Velocity  Pressure  and 
Height  Pressure  for  specific  free-fall  objects  and  dis¬ 
pensers.  In  the  case  of  the  more  complex  ballistic  solu¬ 
tion  (see  Table  1),  the  manuals  include  instead,  the  values 
of  the  Ballistic  Index  for  specific  free-fall  objects.  The 
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manuals  must  also  include  the  Loran  Rectangular  Coordi¬ 
nates  of  the  Loran  transmitters  and  other  basic  chain  in¬ 
formation. 

The  values  of  Time  of  Fall  and  Trail  as  functions  oZ 
Velocity  and  Height  Pressure  are  related  to  the  more  con¬ 
ventional  Time  of  Fall  and  Ballistic  Range  as  functions  of 
True  Air  Speed  and  Altitude,  as  follows: 

At  any  position  in  the  atmosphere, 

HP  is  the  true  static  pressure,  i.  e. ,  Height  Pressure, 

P  is  the  ram  pressure, 
a  is  the  velocity  of  sound, 

TAS  is  the  true  air  speed,  and 
M  is  the  Mach  number. 

Also, 

GP  is  the  static  pressure  at  the  target  point  on  the 
ground. 

The  two  independent  variables  to  be  used  are  the  Ballistic 
Pressure  BP,  defined  as  (GP  -  HP)  and  the  Velocity 
^Pressure  VP,  defined  as  (Pr  -  HP)  (see  sections  on  Air 
Data  and  Air  Release  Point). 

In  isentropic  flow. 


_R=(i+ZJLi.  m2) 

HP  2  ' 


y 

y- 1 


9  3.5 

=  (1  +  0.  2  — ■? — )  for  air. 
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and 


VP  =  HP 


3.5 


(1  +  0.  2 


TAS 


-  1 


These  equations  are  for  use  with  true  Air  Data.  If  the  re¬ 
constituted  ballistic  tables  are  prepared  for  use  with  raw 
"Air  Data, "  i.  e. ,  not  corrected  for  AP  by  calibration  (see 
section  on  "Air  Data"  under  "Aircraft  Systems"),  HP  +  AP 
must  be  substituted  for  HP,  and  VP  -  AP  for  VP  in  the 
above  equations.  If  TF  is  the  time  of  fall  and  BR  the  bal¬ 
listic  range,  then  the  trail,  TR,  is  given  by: 

TR  =  TAS*  TF  -  BR. 


From  the  standard  level-release  tables,  one  may  obtain 
TF  and  BR  in  terms  of  three  parameters:  (a)  true  air 
speed,  (b)  altitude  above  the  target,  and  (c)  target  altitude. 
Height  pressure  and  sound  velocity  as  a  function  of  alti¬ 
tude  may  be  obtained  from  Ref.  2.  Using  the  above  equa¬ 
tions,  the  ballistic  tables  can  be  reconstituted  in  terms  of 
VP  and  HP  (for  various  values  of  GP  if  variations  due  to 
GP  are  significant).  Experience  has  shown  that  ballistic 
data  in  terms  of  VP  and  HP  are  relatively  insensitive  to 
GP.  They  are  also  relatively  insensitive  to  mistakes 
made  in  the  selection  of  the  correct  atmospheric  model 
from  Ref.  2.  It  is  suggested  that  the  atmospheric  model 
used  be  chosen  as  representative  of  the  expected  area  of 
operations. 

The  Technical  Manuals  must  include  pertinent  data 
that  define  the  geometry  of  the  Loran  chains.  These  data 
are  the  group  repetition  period  Pj;  the  coding  delays  Dj  j; 
the  Loran  rectangular  coordinates  of  the  transmitters 
Xj^  j,  Yj  ji  and  the  propagation  velocity  (see  Table  1). 
Pj*and  Dj  j  can  be  taken  directly  from  standard  publica¬ 
tions  on  existing  Loran  chains,  but  Xj?  j,  Yj  j,  and  V  must 
be  determined  from  the  published  long’itudes’and  latitudes 
of  the  transmitters  and  the  associated  coding  delays.  The 
method  for  doing  this  is  given  in  Appendix  C. 
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Planning  Materials 


In  addition  to  standard  flight  planning  materials, 
forms  for  chronologically  organizing  the  ground  computer 
input  are  required.  Each  entry  must  be  addressed  to  en¬ 
able  the  ground  computer  to  process  the  input  data  and 
organize  the  output  consistent  with  Tables  1  and  2.  Prop¬ 
erly  executed  forms  will  be  referred  to  as  the  "Flight 
Program.  " 


Flight  Profile  and  Data 


The  result  of  flight  planning  is  the  definition  of  the 
flight  profile  and  the  associated  data.  Data  recorded  on 
flight  program  forms  described  under  "Planning  Materials" 
constitute  this  definition.  These  data  must  be  consistent 
with  Tables  1  and  2,  except  that  the  horizontal  coordinates 
are  given  in  time  differences;  and  altitude  in  geometric 
height  above  sea  level.  The  ground  computer  converts 
these  inputs  to  the  forms  required  by  Tables  1  and  2.  In 
addition,  the  input  data  must  include  initialization  data 
for  the  ground  computer  subroutines.  These  data  are  the 
sea  level  pressure  and  the  coefficients  of  the  equation  of 
state  of  the  atmosphere.  See  the  section  on  "Height  Pres¬ 
sure"  below. 


Ground  Computer 


The  function  of  the  ground  computer  is  to  convert 
flight  planning  parameters  to  the  required  avionic  parame-’ 
ters.  This  function  is  required  because  the  flight  planning 
parameters  that  are  convenient  to  manipulate  manually 
burden  the  avionic  computer  unduly  if  they  are  used  for 
navigating  the  aircraft.  The  ground  computer  program 
must  recognize  a  library  of  codes  associated  with  the  data 
listed  in  the  Flight  Program  so  that  it  can  properly  pro¬ 
cess  each  type  of  input  according  to  "Data  Type.  "  It  must  ' 
also  include  the  three  subroutines  required  for  "Loran 
Rectangular  Coordinates,"  "Height  Pressure,"  and  "Deriv¬ 
atives.  "  Furthermore,  the  output  format  must  be  con¬ 
sistent  with  the  input  format  of  the  "Aircraft  System.  " 
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Loran  Rectangular  Cor  ~dinates 

i 

The  inputs  to  the  functional  blocks,  "Loran  Rec¬ 
tangular  Coordinates"  of  Fig.  10,  are  the  time  difference 
coordinates  of  waypoints,  target  points,  and  touch-down 
points  as  required  to  specify  the  flight  profile.  The  func¬ 
tion  of  this  biock  is  to  convert  the  time  differences  to 
Xj,  Yj  coordinates  of  Table  2.  This  is  done  by  a  subroutine 
of  theJ ground  computer  program.  This  must  be  done  in  ex¬ 
actly  the  same  way  as  prescribed  for  the  Aircraft  System. 
(See  the  "Loran  Rectangular  Coordinate"  section  under 
"Aircraft  System"  for  the  required  algebra. )  • 

Height  Pressure 

The  inputs  to  the  functional  block  "Height  Press*.- 
are  the  geometric  altitudes  abpve  sea  level  of'waypoints, 
target  points,  or  touch-down  points  as  required.  The  func¬ 
tion  of  "Height  Pressure"  is  to  convert  geometric  altitude 
to  the  HPj  or  GPj  coordinates  of  Table  2  consistent  with 
present  mean  atmospheric  conditions.  This  is  done  as 
follows: 

i. 

HP.  or  GP.  =  P  <1  ->aH)  , 

3  3 

where 

PQ  is  the  current  value  of  atmospheric  pres- 
;  sure  at  sea  level  in  the  area  of  operation, 

a  and  bare  coefficients  Which  define  the  equations 
of  state  of  the  atmosphere,  and 

H  is  the  geometric  height  above  sea  level  of 
waypoints  or  ground  points. 

Sea  level  pressure  and  the  two  coefficients  must  be  en¬ 
tered  into  the  ground  computer  along  with  other  flight  data. 
Their  values  are  obtained  from  "Meteorological  Support. " 


I 
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As  an  example,  if  altitudes  were  given  in  feet,  and 
if  the  latitudes  of  flight  operations  corresponded  to  those 
where  the  ICAO  standard  atmosphere  applies,  the  value 
of  coefficient  a  is  6.  875  x  10“6,  and  of  b  is  5.  256.  A 
study  of  the  standard  atmospheres  of  Ref.  2  indicates  that 
the  values  of  a  and  b  are  not  very  sensitive  to  latitude. 

It  is  expected,  therefore,  that  the  operational  needs  of  the 
ground  system  can  be  fulfilled  by  a  very  limited  number 
of  standard  atmospheres  selected  from  Ref.  2. 

Derivatives 

The  function"Derivatives"  is  performed  if  the  sys¬ 
tem  uses  the  simple  ballistic  solution  referred  to  in  Table 
1  under  subparagraph  c-1.  If  the  more  complex  solution 
(Table  1,  subparagraph  c-2)  is  used,  this  function  is 
omitted. 

The  inputs  to  "Derivatives”  are  values  of  Time  of 
Fall  and  Trail,  paired  with  values  of  Velocity  Pressure 
and  Height  Pressure  selected  from  the  ballistic  tables  of 
the  Technical  Manuals. 

As  before, 

TRq  and  TFq  are  the  values  of  trail  and  time  of 
fall  at  the  planned  reiease, 

and  let 

AVPQ  and  AHPq  be  the  allowable  incremental 

u  changes  in  the  planned  values  of 
VP  and  HP  at  release. 

Then,  select  values  of  TR  and  TF  from  the  ballistic  tables 
corresponding  to  VPg  ±  AVPq  and  HPq  ±  AHPq  and  assign 
index  numbers  as  follows: 
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(VP0  -  AVP0> 

<3 

no 

o 

(VPq  +  A  VP) 

<HPq  -  AHPq) 

— 

TRj,  TFj 

— 

HPo 

TR3‘  TF3 

TR0,  TFq 

TR4,  tf4 

(HP0  +  AHPq) 

— 

tr2.  tf2 

— 

Then 


TRH 


TR2  -  TR 
2AHP 


TFH  = 


TF2  -  TFt 
2AHP 


TRV 


TR^  -  TRg 
2AVP 


TFV  = 


TF  -  TF 
4  3 

2AVP 


The  outputs  of  the  function  "Derivatives"  are  TRg,  TFq, 
TRH,  TFH,  TRV,  and  TFV. 


Flight  Data 

All  flight  data  that  do  not  require  processing  are 
transferred  directly  to  the  medium  by  which  the  Flight 
Input  Data  are  to  be  transferred  to  the  Aircraft  System. 

In  this  instance,  the  ground  computer  serves  only  as  a 
means  of  ensuring  the  correct  organization  of  all  the  data 
to  be  transferred  to  the  Aircraft  System. 
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Flight  Input  Data 

The  flight  input  data  are  the  data  of  Tables  1  and  2. 
For  a  specific  system  these  data  include  only  those  that 
pertain  to  the  options  embodied  in  the  system.  The  actual 
medium  used  for  data  transfer  to  the  Aircraft  System  could 
be  a  recording  such  as  magnetic  tape;  or  if  there  were  a 
requirement  for  the  Ground  System  and  the  Aircraft  Sys¬ 
tem  to  be  interfaced  during  flight,  all  or  part  of  the  trans¬ 
fer  medium  could  be  a  communication  link. 
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5.  SUMMARY 


This  report  has  described  an  accurate,  four¬ 
dimensional  aircraft  navigation  system  in  terms  of  the 
functions  performed  by  various  subsystems.  Mathematical 
descriptions  of  these  functions  are  given  in  consistent 
nomenclature.  Functions  of  subsystems  of  the  supporting 
ground  system  have  been  described  similarly.  In  some 
cases,  the  functional  descriptions  have  included  a  number 
of  options,  one  of  which  can  be  chosen  to  suit  the  need  of 
a  specific  system.  Also,  some  of  the  subsystems  may  not 
be  required  for  some  systems.  For  example,  all  the 
modes  of  operation  shown  in  Fig.  3  are  not  usually  a  re¬ 
quirement  for  a  single  type  of  aircraft.  In  such  a  case 
the  various  modes  to  be  deleted  are  obvious.  Deleted 
aircraft  modes  or  deleted  options  with  respect  to  certain 
aircraft  subsystems  allow  the  deletion  of  their  gr-uind  sys¬ 
tem  counterparts. 
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APPENDIX  A 

RESULTS  OF  STATISTICAL  ANALYSIS  OF 
UDORN  MONITOR  DATA 


The  Southeast  Asia  Loran  chain  is  monitored  at 
Udorn,  Thailand.  At  this  site  two  Loran  timers  are 
operated  in  parallel  from  the  same  antenna.  The  time 
difference  coordinates  of  the  antenna  measured  by  each 
timer  are  recorded  continuously  and  used  to  assess  the 
degree  to  which  the  chain  is  operated  within  tolerances. 
From  time  to  time,  the  U.  S.  Coast  Guard  has  made  these 
records  available  for  analysis  at  A  PL.  The  complete  re¬ 
sults  are  described  in  the  references.  The  mean  and 
standard  deviations,  extracted  from  the  references,  are 
listed  in  the  following  table.  TDX,  TDY,  and  TDZ  refer 
to  time  differences  between  the  master  transmitter  and 
slave  transmiUers  X,  Y,  and  Z,  respectively. 
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APPENDIX  B 

MATHEMATICAL  DEVELOPMENT  OF  THE  NAVIGATION 

COORDINATES 


Let  X.,<P.  -  latitude  and  longitude  of  the  i  station, 

11  i  =  M,  A,  B 

and  Xq,<Pq  -  latitude  and  longitude  of  the  origin, 

and  O.  -  central  angle  from  the  origin  to  i**1  station. 

Now  the  origin  of  the  coordinate  system  Xq,<Pq  is  selected 
such  that  =  ft  ^  =  tl-Q.  The  solution  for  Xq,<Pq  is  as 
follows: 

COS°M  =  COsXM  COS(pM  C0SXO  COSP0 

+  cosXM  sin^M  cosX0  sinpQ  +  sin^  sinXQ,  (1) 

cosO^  =  cosXA  cos<PA  cosXq  C03P0 

+  cosX,  simp.  cosX_  sinP.  +  sinX,  sinX~,  (2) 
A  A  O  O  AO 

cosftB  =  cosXg  cos<PB  cosXQ  cos pQ 

+  cosXB  sirtpg  cosXQ  sino0  +  sinXg  sinXQ.  (?) 

Subtracting  Eq.  (1)  from  Eq.  (2)  and  dividing  by  cosXq, 

O  =  (cosX,  cos<pA  -  cosX,-  coso,,)  cos 

A  A  M  M  O 

+  (cosX^  sirup ^  -  cosXM  sin<PM)  sino0 
+  (sinX,  -  sinX,.)  tanX„  . 

A  M  vJ 

Subtracting  Eq.  (1)  from  Eq.  (3)  and  dividing  by  cosXQ, 

O  =  (cosX,  cosrp  -  cosX,,  C09tP  )  COS0_ 

B  B  M  M  O 

+  (cosX„  siro„  -  cosX,,  sinp„)  sinO.  (5) 


B  B 


M  M 


+  (sinX_  -  sinX,T>  tanX~ 
B  MO 


Preceding  page  blank 
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II 


No\  let 


ai  =  cos*a  cos^A  ‘  COSXM  C0SPM’ 
a2  =  cosXA  sinpA  -  cosXM 

a3  =  sinX^  -  sin^M> 

bl  “  cosXB  COS°B  '  ccsXM  COapM* 
b2  =  COsXB  SiWk  ‘  COsXM  C03PM- 
b3  =  sinXB  ‘  SinXM‘ 

Equations  (4)  and  (5)  can  be  written 


a1  co90Q  +  a2  sinoQ  =  -  ag  tanXQ, 

bl  cos°o  +  bo  sin°o  =  ~  b3  tarV 

Now  multiplying  Eq.  (4’)  by  b0  and  Eq.  (5')  by  a  3  and  sub¬ 
tracting 

(a  b„  -  b  a„)  cosO  +  (a0b0  -  b„a0)  sirx?^  =  O.  (( 


13  13 


2  3  2  3 


From  Eq.  (6)  it  is  clear  that 


..  .  -i  /  V3  ‘  Vi  \ 

°  \a2Vb2a3  / 


Now  from  Eq.  (4')  we  have 


x0  =  tan 


-  -  (a,  COS00 


+  a2  ai“,0))  ' 


Now  that  the  origin  of  our  system  is  defined,  we  shall 
rotate  Xj,c^  into  this  new  coordinate  system.  The 
cartesian  components  of  X^,Oi  are 


r.  =  r 
i  e 


cosX.  COSO.  \ 

1  1  \ 
cosX.  sino.  I 

1  1  / 
sinX.  ' 
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where 


re  is  the  radius  o*  the  earth. 

Now  the  components  of  r!  in  our  new  system  are  given  by 


r! 

1 

where  R^ 


=  R2  R1  ri  * 


is  the  rotation  matrix: 


(10) 


(cosc pQ  si npQ  O 
-siwp0  cos(flQ  O 
O  O  1 


(11a) 


and  R2  is  the  rotation  matrix: 


(cosXQ  O  suiXq 
O  1  O 
-sinXQ  O  cosXq 


) 


(lib) 


Now  denoting  the  components  of  r!  by 


r! 

l 


X! 

l 

Y! 

i 

2! 

l 


(12) 


We  now  define  the  pseudo  latitude  and  pseudo  longi¬ 
tude  of  the  stations  by: 


(13) 


(14) 
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Now  if  the  surface  of  the  earth  is  considered  flat,  then 

e:  =  r:  |(X!  -  XT  +  (<p!  -  <pT  |  ,  (15) 


>f  =r2e  ^(X!  -X')2  +  (<p!  -<P*)2j  , 


where  0^  is  the  arc  length  from  a  point  X',  <p'  to  i*h  station, 
i  =  M,  A,  B.  This  is  equivalent  to  the  small  angle  ap¬ 
proximation 

sinco  =  co  , 

2 

cosoo  =1  -  60  /  2  for  small  co. 

Our  coordinate  transformation  has  assured  us  of  small  co 
(typically  less  than  5°).  For  our  case  co  =  0/ r  . 

Now  letting 

r.  W 

reW 

re  X  =  Y- 

r  v  =  X  . 
e 

We  have  for  i  =  M,  A,  B  from  Eq.  (15), 


4 

>1 

>1 


Equations  (17)  and  (18)  can  be  written  as 

<6A  '  eM  +  V"  *  <XA  -  X)2  +  <YA  -  Y)2  ’  (1,,) 


=  (XM  -  X)2  +  (YM  -  Y)2  , 

(16) 

=  (XA  -  x)2  +  (Ya  -  Y)2  , 

(17) 

=  (Xfi  -  X)2  +  (Yb  -  Y)2  . 

(18) 
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<»B  '  8M  +  6M)2  =  (XB  -  X)2  +  <YB  •  Y'2 


(18') 


Now  let 


**A  =  ®A  '  9M  * 

=  SB  "  ®M  * 

Note  that  our  measurements  are  linear  functions  of  Ja 
and  uB  so  uA  and  uB  are  inputs.  Equations  (17  )  and 
(18')  are 

<4  +  2ma  6m  +  6m  =  xa  ‘  2X  xa+x2 


+  <-2yya  +  y2' 


4+2^AeM  +  8M=4-2XXB+X" 


‘A  M  •  M  B 


2  ,  .2 

+  Yi-2YYB  +  Y  • 


Now  subtracting  Eq.  (16)tromEqs.  (19)  and  (20), 


4+-2“A6MsXi+YA+2(XM-XA,X 

+  2<YM  '  V  Y  ‘  <XM  +  Ym' ’  ' 

4+2**BeM=x2B+YB  +  2(XM-XB)X 


+  2(YM  *  V  Y  '  <XM  +  YM>  ' 
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Now  let 

a  =  ^  <1/2  D|»2  -  1/2  C«|)  , 

8=5(dMa-c»b>- 

y.i<l/2A«|-l/2Bu^), 

«  =Z<a“b  -  • 

Equations  (23’)  and  (24')  can  be  written  as 

y  =  «  +  00m* 


X=V+69m. 


(25) 

(26) 


Now  substituting  X,  Y  from  Eqs.  (25)  and  (26)  into  Eq.  (16) 
gives 


Kn  s  (Xa/t 

M  M 


1  -  66m)2  +<*„-«-  ^M)2  • 


M' 


M 


(27) 


Expanding  Eq.  (27), 


9m  =  ym  +  4-2Ym“-2Xm’' 


'M 

Collecting  terms 

,2  .  »  2  ...2 


-2YM^M-2XM  66M+“2+y2 

+  20fse„+/s2e2 +2yee  +62e^. 


o  =  oz  +  62  - 1)  ej,  +  2[<wi  -  vM/[  +r6  -  xM  6)  eM 


M 

+  (YM-rf2  +  (XM-y)2. 
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APPENDIX  C 

DETERMINATION  OF  TRANSMITTER  LOCATIONS  IN 
LORAN  RECTANGULAR  COORDINATES 


Appendix  B  develops  the  coordinate  transformation 
that  is  to  be  used  for  aircraft  navigation.  For  simplifica¬ 
tion  of  the  mathematics,  the  development  is  based  on  a 
spherical  model  of  the  earth.  The  resulting  transforma¬ 
tion,  referred  to  as  "Loran  Rectangular  Coordinates,  " 
adequately  fulfills  the  requirements  for  accurate  naviga¬ 
tion  as  described  in  the  text  of  the  report.  However,  com¬ 
plications  develop,  principally  in  the  vicinity  of  baseline 
extensions,  if  the  Loran  transmitters  are  not  located  in 
Loran  Rectangular  Coordinates  with  sufficient  accuracy 
Analysis  of  this  problem  has  shown  that  a  spherical  earth 
model  is  not  sufficient  for  locating  the  transmitters. 

Specifically,  the  problem  is  that  the  locations  of 
the  transmitters  that  would  be  obtained  by  the  procedures 
of  Appendix  B  do  not  constrain  the  Loran  Rectangular 
Coordinate  transformation  to  exist  everywhere  for  '  real- 
world"  measured  values  of  time  differences.  Tl.  mathe¬ 
matics  allow  imaginary  values  of  Xj*  and  Yr  under  some 
circumstances.  The  desirable  goal  of  precluding  these 
circumstances  may  be  achieved  using  "corrected"  coor¬ 
dinates  of  the  transmitters  and  propagation  velocity  while 
navigating.  The  correct  values  may  be  obtained  by  per¬ 
turbing  the  transmitter  locations  and  the  ve’ocity  of  propa¬ 
gation  in  such  a  way  as  to  match  the  time  differences  that 
identify  the  six  baseline  extensions  to  the  real  values. 

The  procedure  for  doing  this  is  as  follows: 

1.  Obtain  published  values  of  the  geodetic  loca¬ 
tions  of  the  transmitters  and  compute  Loran 
Rectangular  Coordinates  as  Appendix  B. 
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2.  Obtain  the  real  baseline  lengths  in  microseconds 
from  calibration  data. 


Then,  using  the  transmitter  locations  and  baseline 
lengths,  set  ud  the  live  constraining  equations  in  the  seven 
unknowns.  These  unknowns  are  the  amounts  by  which  the 
three  pairs  of  transmitter  coordinates  and  the  propagation 
velocity  must  be  perturbed.  Solve  these  equations  for  five 
of  the  unknowns  in  terms  of  the  constants  and  the  other  two 
unknowns.  Obtain  a  minimum  variance  solution  for  these 
two  unknowns  and,  using  these  values,  obtain  the  values 
of  the  other  five  unknowns.  The  constraining  equations  are 

as  follows: 

(X.  -  X.  +  AX.  -  AX.)2  +  (Y.  -  Y.  +  AY.  -  AY.)2 

=  (V  +  AV)2  (L.  .)2  , 

*■»  3 

where 

i  =  ?/.,  A,  B, 
j  =  M,  A,  B, 

i  ^  3  in  the  same  equation, 

X,  Y  are  Loran  Rectangular  Coordinates 
V  is  the  propagation  velocity,  and 
L  r  re  the  baseline  lengths  in  microseconds. 

i;3 

Permutations  of  the  indices  in  the  equation  above  result 
in  three  equations.  Also, 

(X.  +  AX.)2  +  (Y  +  AY  )2 

i  i  i  1 

=  (X.  +  AX.)2 +(Y  +AY)2. 

3  3  3  3 
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Permutations  of  the  indices  result  in  three  of  these  equa¬ 
tions,  but  one  of  the  three  is  redundant.  Also, 


SAX  2  +  SAY.2  +  AV2  =  minimum, 
i  1  v 


where  R.  =  Vx;2  +  using  the  values  obtained  in 
Appendix  B.  This  is  the  function  to  be  minimized  by  a 
least  square  computational  routine. 


These  equations  can  be  most  easily  solved  by 
linearizing  where  appropriate  and  iterating  th®  soluUom 
Starting  values  are  obtained  by  Appendix  B.  Final  values 
are  ^"corrected”  values  to  be  used  in  the  navigation 

equations. 
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